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Abstract 

We review a surprising correspondence between certain two-dimensional integrable 
models and the spectral theory of ordinary differential equations. Particular em- 
phasis is given to the relevance of this correspondence to certain problems in VT- 
symmetric quantum mechanics. 
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1 Prelude 



This short review is about a surprising hnk between integrable models and ordinary 
differential equations |21 El 13 ' relevance in the study of the novel kind of non- 
Hermitian quantum mechanics studied by Carl Bender and many others [SlIHl 171 IHIIUI ITUl 
El Iia Iia ini lEl lEl ini CHI lEl EOl ini 123 application of techniques from the 

theory of integrable models to answer a long-standing question concerning the reality 
of the energy levels of a particular set of non-Hermitian operators. 

In the theory of ordinary differential equations the adjective "exactly solvable" usu- 
ally indicates that the problem is fully soluble, with solutions expressible in terms of 
some previously-defined mathematical functions. In the context of spectral theory of 
second-order differential equations the classical example is a Schrodinger problem with 
harmonic potential: 

- tp'^x) + = En^{x) , xp{x) G L^{M) . (1.1) 

This equation is mathematically soluble in the sense that the eigenfunctions are given in 
terms of elementary functions. On the other hand, a mathematician would be reluctant 
to call a Schrodinger equation with potential with a G IR"^ exactly solvable. However 
we shall shortly see that this family of quantum systems is somehow equivalent to a two 
dimensional exactly solvable model: the six-vertex model in its scaling limit. There is 
clearly a mismatch between the definitions of exact solvability in the two fields. Naively 
speaking, in two dimensions there is a special class of statistical models for which at least 
one physically-interesting quantity, the density of free-energy in the thermodynamic 
limit, can be computed exactly. For this reason these class of systems are given the 
name of "exactly solvable" or "integrable" models. (The existence of a complete set 
of independent, commuting conservation laws is a more precise characterisation of an 
integrable system, but a full discussion of this would go well beyond the scope of this 
review.) 

Notice that this does not in general mean that other quantities, such as, for example, 
the finite lattice free energy or the correlation functions, can be exactly determined. This 
sounds like a negative note (and it is!!), but during more than 30 years of integrability 
the scientific community has partially overcome this negative aspect by introducing a 
number of powerful tools which allow an integrable system to be studied in great detail 
- works particularly relevant to the current story include [231 1241 125j . This review is 
about the use of just one these tools in the ODE context. Of course many techniques 
have also been developed in the ODE framework (see for example |261 1271 128j). and 
integrable models can, and to some extent already do [21123) profit from these as well. 

2 The IM side: the six-vertex model 

Let us start from the integrable model side, and consider an x M lattice model with 
periodic boundary conditions and, for an irrelevant technical reason, N/2 even. On each 
link of the lattice we place a spin taking one of two values: this is conveniently denoted 
by placing arrows on the links, as in figure ^ 
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Figure n A configuration in the six- vertex 
model. 
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Figure 12 The t'i's are real for the ground 
state. 



We shall constrain the set of configurations as follows: 



• Only those configurations of spins which preserve the 'flux' of arrows through each 
vertex are allowed; 

• We only consider the 'zero field' six-vertex model, which has an additional '4-spin 
reversal' symmetry implying the model is invariant under simultaneous reversal 
of all arrows. 

Thus the 'local Boltzmann weights', numbers assigned to each vertex depending on the 
spins next to that vertex, can be parameterised in terms of just three quantities: 
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In fact the overall normalisation factors out trivially from all calculations, and the 
remaining two degrees of freedom can be parametrised using the two variables 



the spectral parameter , 77 : the anisotropy 



(2.1) 



as 



a = sinh(ir/ — v) , h = sinh(ir/ + v) , c = sinh{2iri) . (2-2) 

The relative probability of finding any given configuration is simply given by the product 
of the Boltzmann weights at each vertex. The six-vertex model introduced above is a 
classic example of an integrable lattice model. 

The first quantity of interest - the partition function Z - is the sum of these numbers 
over all possible configurations: 



(2.3) 
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One popular technique to calculate Z makes use of the so-called transfer matrix T, 
which performs the sum over one set of horizontal links: 
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(2.4) 



In terms of T the partition function is given by 

Z = Trace [T^] . 

In the limit M ^ oo with finite the free energy per site can then be obtained as 

/ = j4 In Z = j4 In Trace [T*^] ~ ^ In to , (2.5) 

where to is the ground-state eigenvalue of T. The problem is thus reduced to the 
determination of to ; and since the model is integrable, there are many methods to 
achieve this end. The detailed description of these methods goes beyond the scope of 
this review, but it is worth mentioning that they usually lead to a set of non-linear 
constraints (functional relations or non-linear integral equations) on to(z^), or, more 
simply, constraints on the zeroes of a related function qQ{v)., which we shall introduce 
shortly. These constraints are very powerful even at finite N but they lead to a closed 
expression for the free energy only in the thermodynamic limit iVf, N ^ oo. 

Here we shall sketch the logical flow of a method developed by Baxter (for more 
details, see j23j ) leading to a relation known as the TQ-system that can also be easily 
deduced starting from the ODE side. This should provide a first clear hint of what we 
mean by an ODE/IM correspondence. 

Baxter began by showing that there exists an auxiliary function qo{i^), 



n-l 



sinh(zy - i^i) , 



(2.6) 



such that the following 'TQ relation' holds 

to{i^)qo{^) = a{u, ri)^qQ{v + 2ir]) + h{v, rj)^ qo{i^ - 2ir]) 



(2.7) 



This is perhaps a puzzling step to take: we want to find to(z^) and we now claim that the 
relation (|2.7|) . which just defines to(z^) in terms of another unknown function qoii^), will 
somehow help us in this task. The explanation is simple: the constraint ()2.7() should be 
combined with a knowledge of the analytic properties of both tg and qo. In particular 
the fact that to(j^) and ^o(z^) are entire functions of means that a powerful constraint, 
the Bethe ansatz system, is immediate. By definition 



and combining this with the relation (|2.7)) and the fact that to(z^) is entire leads to the 
following set of Bethe ansatz equations (BAE): 



a^{ui,ri) qo{ui + 2ir]) 
b^it^i,v) Qoii^i - 2ir]) 



0,l,...,n-l . 



(2S 



This already looks to be a more serious constraint than 1)2. 7|) . but this is not yet sufficient 
since (|2.9j) has many sets of solutions {fj}. However, one can argue that there is a 
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particular solution, with n = N/2 and all of the Ui real (depicted in figure which 
uniquely corresponds to the ground-state eigenvalue of T, i.e. Iq. 

To prepare the ground for the connection with differential equations, let us define 



Aj = e""' , uj = —e 
so that the Bethe ansatz equations 1)2. 9|) become 



i2il 
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n=0 



(A„ - XiUJ 2) 



(2.10) 



(2.11) 



We should mention one final detail: currently correspondences with ordinary differen- 
tial equations have only been established with integrable models in their continuum, 
conformal limits. For the six-vertex model this limit is found by sending 







N —>■ oo and a = e" 
with aN kept finite. Sending 

a ^ as a ^ a/N 
one discovers that the Aj's with i <^lnN rescale to zero as 

A,- ~ K-a^''/^ ~ . 



(2.12) 



(2.13) 



(2.14) 



For 7r/4 < rj < Tr/2 the limiting product converges with no need for extra regulating 
factors, and one obtains 
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n 

n=0 



v„ — EiUJ 

{En — EiUJ~ 



(2.15) 



Finally, the six-vertex model can be generalised to incorporate twisted boundary con- 
ditions without losing integrability. The twist is introducing by making the following 
replacement in the definition of the original model: 
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(2.16) 



In the conformal limit described above the introduction of the twist leads to the more 
general BAE 
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n 

n=0 



{En - EjLo') 
{En - EiOJ-^) 



(2.17) 



We shall shortly see how to derive the same set of equations starting from a differential 
equation. 



3 T^T-symmetric quantum mechanics 

Before we go into the details of the precise link between ODEs and integrable models 
we shall first give a description of the non-Hermitian quantum mechanical problems 
mentioned in the prelude. 
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Researchers working on integrable 1+1 dimensional massive quantum field theories 
have become quite accustomed to systems described by non-Hermitian Hamiltonians 
which nevertheless possess, at least in a certain range of the parameters, a real and 
positive energy spectrum. A standard example, which has been studied since the mid- 
1990s, is the scaling Lee- Yang model (see, for example, [201 )• This is the theory of 
a massive relativistic bosonic field (/> with cubic interaction gcjr". Along with the A(/)^ 
case, it is one of the standard text-book examples used for practising Feynman diagram 
technology [HJ. However, since the potential is unbounded from below, the Lee- Yang 
model with real g suffers from obvious pathologies. Surprisingly, at least in the 1+1 
dimensional integrable case, these pathologies disappear in the fully renormalised scaling 
limit at g = i\g\. We invite the interested reader to consult |3U1 1321 15^ to appreciate 
various aspects of this phenomenon, and to note in particular the striking similarity 
between the plots in figures 6 and 7 in Ref. |SS1 and figure IHK below, which shows the 
energy spectrum of a non-Hermitian quantum-mechanical problem, (|3.2|) . 

Indeed, it was the study of the 1+1 dimensional Lee- Yang model that led Bessis and 
Zinn-Justin to the following, rather simpler, question: is the energy spectrum associated 
with the Schrodinger equation 

(f 

- T^i^ix) + ix^i^{x) = Enip{x) , ^{x) E L2(]R) (3.1) 

real? Perturbative and numerical studies led Bessis and Zinn-Justin to conjecture that 
the spectrum {En} is indeed real, and positive But how does one prove this 

statement analytically? 

Later, in 1997, Bender and Boettcher [7j considered the spectrum of the following 
generalisation of Bessis and Zinn-Justin's problem: 

- -T^i^ix) - (ix)^ij(x) = EnMx) , (N real, > 0) (3.2) 

with the associated boundary conditions to be specified shortly. Again, extensive nu- 
merical checks led the authors of to a precise reality conjecture. They also remarked 
that both equations p.lj) and (|3.2|1 are invariant under the action of the operator VT: 

V (parity reflection): x —x , p —p , 

T (time reversal): x ^ x , p — > —p , i ^ —i . 

The special class of non-hermitian quantum systems with VT symmetry has only re- 
cently started to be studied in depth. These theories are mathematically interesting, 
and they may have an important role in physics I35J . 

Returning to the ODE (|3.2|) . we note that there appear to be a couple of problems 
with the generalisation. Firstly, for non-integer values of N the 'potential' — (zx)^ is 
not single- valued and we need to add a branch cut, which we shall put along the positive 
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imaginary x-axis: 



Im 



Re 



Figure 13 The branch cut, and a ray in the complex x-plane. 

Secondly, when N reaches 4 the 'potential' is — x^, and the eigenvalue problem runs 
into difficulties since all the solutions to ()3.2|1 decay algebraically as |x| — > oo on the 
real axis. To overcome this problem it is necessary to enlarge the perspective and treat 
x as a genuinely complex variable. Consider solutions to the ODE in the complex x- 
plane. For most values of arg x , up to an overall normalisation there will be a unique 
exponentially-decaying solution at large \x\ (this solution is called "subdominant" ) and 
any other linearly independent solution will be exponentially growing (these solutions 
are called "dominant"). However, whenever 

vr Svr Svr 

all solutions decay algebraically. The lines along which this occurs are known as 'anti- 
Stokes lines' and the wedges between the anti-Stokes lines are the 'Stokes sectors'. Since 
all functions involved are analytic, we can continue the wavefunction along some other 
contour in the complex plane, and as long as it does not cross any anti-Stokes lines 
the spectrum will remain the same. At = 4 an anti-Stokes line coincides with the 
real axis and the correct analytic continuation of the original problem is achieved by 
bending the wavefunction contour down into the complex plane, as shown in figure 0] 

We conclude from this discussion that each pair of Stokes sectors in which the 
wavefunction is required to decay at large \x\ defines a different eigenvalue problem. In 
the terminology of the WKB method, these are related to 'lateral' connection problems. 
Contours which instead join x = to x = oo lead to what are called 'radial' (or 
'central') connection problems. Figure El depicts a sample of the possible lateral and 
radial wavefunction contours. 




Figure 01 A possible wavefunction Figure [31 Some further quantization 
contour for > 4. contours. 



Questions in "PT-symmetric quantum mechanics are all related to lateral problems, 
with one particular pair of Stokes sectors selected. The first eigenproblem to arise in 



6 



the ODE/IM correspondence was of radial type with the contour defined along the 
positive real axis. However it turns out that the PT-symmetric quantum mechanical 
problems (||-i.2|) are intimately related with certain radial problems, in a way that will 
be discussed in the next-but-one section. 

4 Numerical evidence 

We now present some numerical evidence concerning PT-symmetric problems of 1)3. 2|) . 
In figure IHK part of the spectrum of (|3.2|) is plotted. This plot, together with figures 
inb-d, was obtained in jlUj via a non-linear integral equation. The results shown in the 
first plot had previously been obtained in [71 |S] by a direct numerical treatment of the 
differential equation in the complex plane. 




0,5 1 1.5 2 0.5 1 1.5 2 

M M 



Fig.Et: / = -0.001. Fig. Eli: I = 0.001. 

FigureEl Eigenvalues of the Hamiltonian — (ia;)^*^ + l{l+l)x^'^ . 

Prom figure EKi for N = 2M > 2 the spectrum is, within our numerical precision, 
real and positive, while as moves below 2 an infinite number of energy levels pair 
off and become complex. The transition to infinitely-many complex eigenvalues was 
interpreted in Pj as a spontaneous breaking oiVT symmetry*. 

'Since \VT,H] = it follow that if an eigenfunction ^l/ of // is also an eigenfunction of VT with 
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We shall tackle the reality question for N > 2 analytically. For this purpose it 
is convenient to enlarge the perspective by including two extra parameters, a and I, 
a generalisation that does not add any extra technical difficulties but gives a wider 
phenomenology. It is also convenient to trade the parameter A'^ for M = N/2. The 
more general theory is |lfl| lllj 

- ^V'(^) - {{ixf'' + a{ix)''~^ + fcll)V'(:E) = En^P{x) . (4.1) 

Even with a = 0, for — 1 < Z < the additional angular-momentum term has a remark- 
able effect on the connectivity of the spectrum, as can be seen in the middle two plots of 
figure ini Moreover, we shall show in the next sections, as in |111I12) . that the spectrum 
of (gU) is 

• real if a < M + 1 + \2l+l\ 

• positive if a < M + I - \2l+l\ . 

The proof makes use of integrable model technology. While there will not be space to 
go into details below, we remark that these ideas can also be used to study the way that 
the energy levels merge to become complex LSHj . 



5 Analysis of the Schrodinger equation 

As in |in| lllj. we start the analysis by considering a related differential equation 

with X and E possibly complex and M > 0. 

We shall need a couple of important facts. Firstly, the equation has a solution 
y = y{x, a, E, I) such that: 



y is entire in E and x (though, due to the branch point at x = 0, x must in general 
be considered to live on a suitable cover of the punctured complex plane) ^ 



for M > 1, y and y' = dy/dx admit the following asymptotic representations 



-M/2-a/2 M/2~aj2 

y- ^ eM-jki^^'^') , y'~-^^=^exp(-^x-+^) (5.2) 



for |x| — > oo in any closed sector contained in the sector | argx | < 2MT2 (though 
extra terms appear for < Af < 1.) 



eigenvalue A then PT i?* = _B*A* = H PT* = BA* and E = E* . 

^The entirety of y was first proved by Sibuya (see Ref. |26|V His work concerned only the case 
I = a = Q, 2M £ N, but the result also holds for the more general situation of eq. (15. It . so long as 
the branching at the origin is taken into account. In this respect we should mention that the I = 0, 
a = 0, 2M £ case was explicitly discussed by Tabara in |37) . while the generalisation to a potential 
P{x)/x^ with P{x) a polynomial in x was studied by MuUin |38|. and more recently in |39| . It is also 
worth noting that with a change of variable it is possible to map eq. (15.111 with a £ TR, I £ TR, 2M £ Q"*" 
onto particular cases of those treated in |38| . 
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Secondly, setting x = pe*^ as in figure|31with p real, and denoting the sector 
'^"^ figure [71 we see that, as p ^ 00 with 9 fixed, 

• in 5o : y ^ (y is subdominant in Sq) 



at 



■ 2M+2 



y decays algebraically 



in S±i : y ^ 00 (y is dominant in S±i). 



k-K 



M+1 



< 



S_3= S 




Figure [71 The sectors for the potential x^. 

It is then easy to see that the uniquely-determined solution y{x, E, a, I) defines an 
associated set of functions 



Vk 



(5.3) 



which are also solutions of H5.1|) for integer A;, and any pair {yu^yk+i} forms a basis of 
solutions. We can therefore write y_i as a linear combination of yo ^-nd yi. The result 
has the form 



T{E, a, l)yo{x, E, a, I) = y-i{x, E, a, I) + yi{x, E, a, I) , 



(5.4) 



where the function T is called a Stokes multiplier. Keeping 5Re/ > —1/2, the leading 
behaviour of y near x = at generic E is 



y{x,E,a,l) Q{E,a,l)x ' + . . . , 
and in terms of the shorthand notation 

= T{E, ±a, I) , Q W = Q W {E) = Q{E, ±a, /) 



(5.5) 



(5.6) 



we find the following relations which intertwine the as-yet undetermined functions T^^^ 
and QW: 

T^+\e)Q^+\e) = uj-^^^^ Q^-\u^^'^ E) + uj^^^^ Q^-\uj-^^'^ E) (5.7) 
T^-\e)Q^-\e) = uj-^^^^ Q^+\oj^^^ E) + oj^^^^ Q^+\oj-^^'' E) . (5.8) 

Notice the striking similarity between equations (|5.7|) . (|5.8jl and the Baxter TQ-system 
of equation H2.7|) . 

At the zeroes {£'fc(/,Q)} of Q{E), the leading behaviour of y at the origin changes 

to 

y(x, ^fe, a, ~ Q{Ek,a, -l-l)x^+^ + ... , (5.9) 



and y{x,E,a,l) decays at the origin as well as at infinity. This implies that Q{E) is 
the spectral determinant encoding the eigenvalues of 1)5. If) for boundary conditions of 
radial type. For M > 1, we can use the Hadamard factorisation theorem to write Q as: 

Q{E,l,a)=Q{^,l,a)\{[l- — \ . (5.10) 

n=0 

Both r(±)(E) and QW(-E) are entire in E, so the LHS of the relevant TQ equation 
((|n3) or (|^ ) vanishes at 

E = E^^^ = Ek{l,±a) , (5.11) 

and the following system of equations of Bethe ansatz type for the energy spectrum is 
obtained 

°° / , ,-2M \ 

_ , -21-1-cy. 




En 
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— U)~ 


-2M 


Ek 
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— U) 


2M J^i+) 




) 


— UJ~ 


-2M Tpi-) 


En 





— UJ 


2M 



uj-''-'+'^ . (5.12) 

n=0 \ En ' - UJ '^^^^ E^, ' ) 

At a = 0, {e'^^^^ = {E^ ^} and the set of equations ()5.12|) reduces to the Bethe ansatz 
system (|2.17p of the six-vertex model in its continuum limit. (For a ^ 0, the mapping 
is instead to a theory called the three-state Perk-Schulz model pil40j.^ 
The Stokes multiplier 



y-iix) yi{x) 
y'-i{x) y[{x) 



(5.13) 



T{E,a,l) = W[y-i,yi] = Det 

given here Wronskian, vanishes if and only if 

T^[y_i,yi] =0 y-i & yi are linearly dependent. (5-14) 

This holds if and only if the ODE has a solution decaying in the two sectors S^i 
and Si simultaneously. Since (|5.1j) is related to the PT-symmetric problem (|4.ip by 
the transformation x x/i, E — > —E, this means that T{—E, —a, I) is precisely the 
spectral determinant for the generalised Bender-Boettcher problem (|4.1j) . Thus the TQ 
relations (|5.7|) . (|5.8|) encode the spectra of both radial and lateral eigenproblems. 

6 A simple proof of the reality property 

In order to prove the reality and positivity claims made at the end of section@]we return 
to equation (|5.7|) 

and define the zeroes of T^^^ = T{E, a, I) to be G {— Afc}. Setting E = —\k and using 
the factorised form for Q*-"-* we get 

°° / _l_ , ,-2M X \ 

n V) ,M = -^"'""^ ' ^ = 0' 1' ■ ■ • ■ (6-2) 
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Since the original eigenproblem is invariant under I —l—l we can assume I > 
— 1/2. Then each En ^ is an eigenvalue of a Hermitian eigenproblem associated with 
TC{M, —a,l), and hence is real. It is also easy to show that the eigenvalues are all 
positive, provided a < M+2/+2 Taking the modulus^, using the reality of the 



\ and writing 

oo 

ni 

n=0 



= |Afc| exp(z(5fc), we have 
■(£;i-))2 + |Afc|2 + 2Ei-)|Afe|cos( 



27r 
A/+1 



|Afc|2 + 2Ei ^|Afc|cos( 



2-K 
A/+1 



1. 



(6.3) 



For a < M + 2/ + 2, all of the E\i are strictly positive and each single term in the 
product is either greater than, smaller than or equal to one depending only on the cosine 
terms. To match the RHS we therefore must have 



cos 



2-K 

M+1 

Since M > 1 the only possibility is 



cos( 



2-K 
Af+1 



nvr . 



n G Z 



(6.4) 



(6.5) 



and the eigenvalues are indeed real. Relaxing the condition on / we have shown the 
spectrum of 1)4. 1|) is 



real 



if a<M + l + |2/+l| 



and using continuity in M to keep track of the signs of the eigenvalues (see |11| ) it can 
then be seen that the spectrum is 

• positive if a < M + 1 - |2/+1| . 

Referring to figure |B1 the spectrum is entirely real for (a, /) G BUCUD and positive for 
{a, I) G D . For M < 1, the order of Q^~^ is greater than one, the Hadamard-factorised 
form for Q^~\E) no longer has such a simple form, and the proof breaks down. In fact 
we know that most of the A^'s become complex in this region. 

Finally, we remark that while the above constraints on the parameters M, a and I 
are sufficient they are not necessary, as can be seen by studying figure |H1 for the case 
of M = 3. The full domain of unreality obtained numerically in J2] is shown as the 
interior of the curved line, a proper subset of A. 
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Figure IHl The 'phase diagram' at 
fixed M. 




Figure El The domain of unreality A 
for M — 3 obtained numerically. 
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7 Conclusions 



We hope to have shown that there is a very interesting relationship between the con- 
formal limit of two dimensional integrable models and the spectral theory of ordinary 
differential equations. The instance described in this review involves the six-vertex 
model in its conformal limit and an interesting class of PT-symmetric quantum me- 
chanical systems with complex potential. As an application of this correspondence we 
have briefly sketched the proof of a conjecture due to Bessis, Zinn-Justin, Bender and 
Boettcher concerning the reality of the spectrum of a particular class of PT-symmetric 
operators. 
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